Statistics of tumbling of a single polymer molecule in shear flow. 
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We present experimental results on statistics of polymer orientation angles relatively to shear 
plane and tumbling times in shear flow with thermal noise. Strong deviation of probability distri- 
bution functions (PDF) of these parameters from Gaussian was observed and a good accord with 
theory was found. The scaling relations of PDF widths for both angles as a function of the control 
parameter Wi are verified and compared with numerics. An universal exponential PDF tail for the 
tumbling times and its predicted scaling with Wi are also tested experimentally against numerics. 
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Dynamics and statistics of a single polymer molecule 
in stationary as well as in random flows have recently 
attracted attention of both experimentalists [1-4] and 
theorists [5-10]. Stretching dynamics and statistics and 
coil-stretch transition in these flows were investigated in 
detail. Besides, another remarkable effect was first ob- 
served in a shear flow [2], namely large fluctuations in a 
polymer elongation due to end-to-end aperiodic tumbling 
(see Fig. ID). 

Recently rather extensive theoretical [11,12] and nu- 
merical [13,14] efforts were conducted with the goal to 
understand the statistics of the angular orientation of 
a polymer molecule and of the tumbling time in a ran- 
dom velocity field with a mean shear. Shear flow with a 
thermal noise was considered there as a particular case. 
A role of thermal noise on a solid rod tumbling in a 
shear flow was first considered in Ref. [15]. Then the 
role of the Brownian fluctuations in the polymer dynam- 
ics and statistics was studied in numerical simulations 
[7], where, however, only the power spectral density and 
the statistics of polymer extension in a shear flow were 
investigated. An average polymer extension and angu- 
lar orientation were theoretically considered also in Ref. 
[8] , where the results of calculations were compared with 
light scattering measurements [16,17]. 

In this Letter we concentrate on statistics of angular 
orientation and tumbling and scaling relations of their 
characteristics for a single DNA molecule in a shear flow, 
when thermal fluctuations are the main cause for tum- 
bling. 

It is well known that a probability distribution func- 
tion (PDF) of the end-to-end vector R for a polymer de- 
scribed by a dumbbell model with a linear relaxation in 
a simple shear is Gaussian [7]. Nevertheless, PDFs of the 
polymer extension, \R\ = R, as well as polymer angular 
orientation are strongly non-Gaussian due to anisotropy 
introduced by the shear. A functional form of the PDF of 
the polymer extension in a shear flow was first identified 
experimentally [2] and then explained theoretically and 
numerically [7]. 

The main result of the recent theory [11,12] is the pre- 



diction of the intermittent (non-Gaussian) statistics of 
polymer angular orientation in a shear flow. Moreover, 
it was also shown that intermediate asymptotic of the 
angular polymer statistics bears some universal features 
independent of the nature of polymer random excitation 
[11,13,14]. 

The universal features of the orientation statistics are 
related to a deterministic nature of the angle evolution at 
large angles, where a random pumping can be neglected. 
Then the theory [11,12] predicts the following results on 
the angular distributions and tumbling time statistics: 
(i) at \4>\ > A(j) ~ WT 1 one gets P(<j>) oc (sin (/)) , 
where A4> is the width of P(4>); (ii) at 1 > \0\ > A0 the 
PDF tail for 9 is determined as P(6) oc |#|~ 2 ; both scal- 
ings follow from the deterministic process; (iii) the tail of 
PDF of tumbling time between two consequent flips, r^, 
is exponential at ^> r t , where r t is the characteristic 
tumbling time that is proportional to r re i . Thus, the ma- 
jor theoretical finding is the intermittency phenomenon 
in the statistics of a single polymer tumbling that results 
in either algebraic or exponential tails of the distribution 
of polymer angular orientations and tumbling times, re- 
spectively. 

Another theoretical result verified by extensive numeri- 
cal simulations based on linear dumbbell as well as FENE 
models describes the scaling relations for the width of 
the angular distributions and the characteristic tumbling 
time with Wi [11-14]. The scaling can be derived by sim- 
ple physical arguments used by the theory [11,14]. The 
dynamic equation for one of the molecule orientation an- 
gles (see Fig.lA), namely <fi, in the region of <fi <C 1 can be 
written as j& = -s sin 2 Or e z)~ 1/2 ?70 [11-14], where 

R g is the polymer gyration radius, and is the white 
thermal noise. In a stationary state one gets from a bal- 
ance of shear and noise terms the following estimates for 
the rms of the angular fluctuations (or the PDF width), 
A</>, and the characteristic tumbling time, r t : 



A<t>ocWi-i& 
R 



r t oc T rel Wi '(%■) 



(1) 



Then one can consider two limits at Wi » 1 (or small 
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A<f> and AO). In the linear extension limit, R -C Rmax: 
where R ma x is the maximum polymer extension, a typi- 
cal polymer elongation, R, can be estimated as R g ■ Wi 
that is realized at sufficiently small extensions (up to 20% 
of R m ax [1,4]). It leads to [11-14] 

A(j) oc Wi' 1 , T t <XT re i. (2) 

In the opposite limit of a non-linear regime, where typi- 
cally Wi > Rmax/Rg, a stretched polymer can be treated 
as a rigid rod with R ~ R m ax subjected to the Brownian 
motion [15]. Then the resulting scaling appears to be 
strikingly different [11-14] : 

A<j> oc wr 1/3 , T t oc T re iWr 2/3 . (3) 

Similar scalings can be also expected for AO. 

The experiments were carried out in a shear flow in 
two flow configurations. To measure we used a shear 
flow in a gap between the flat bottom of the uniformly 
rotating glass rod of radius r\ — 1.5 mm and a cover 
slip (see for details of the set-up [4]). The measurements 
were made at the fixed height of ~ 30/xm above the cover 
glass in the region located at radius r ~ 700/mi, where 
the main shear occurs in the vertical plane. Then only 
can be measured in this configuration (see Fig. IB). 

To measure <f> a shear boundary layer of a Poiseuillc 
flow in a micro-channel was used. The micro-channel 
of ~ 500/mi wide and ~ 100/mi deep was produced 
from a single cast of a silicon elastomer Polydimethyl- 
siloxane (PDMS) using the soft lithography method [18]. 
The measurements were carried out at about the middle 
height of the channel in a small region at a distance of 
~ 30/im from the channel wall, where only <f> was detected 
via the cover glass (see Fig.lC). 

Particle image velocimetry (PIV) measurements were 
conducted with 0.2/mi fluorescent beads in both flow con- 
figurations. In the swirling flow the shear rate in the ver- 
tical plane is rather constant at the location of measure- 
ments of ~ 30/im (Fig.2A), while the azimuthal veloc- 
ity in the horizontal plane almost independent of radius 
in the region of measurements (inset in Fig.2A). In the 
channel flow the longitudinal velocity measured in the 
horizontal plane changes rather linear in the observation 
region of ~ 30/im from the wall, while looks fairly con- 
stant in the vertical plane at the middle height of the 
channel within the focus depth (0.4 /mi) (see Fig.2B). 
In both cases the measured velocity profiles are com- 
pared with numerical calculations, and good quantitative 
agreement is found (solid lines, Fig.2A,B). A limitation 
to conduct PIV measurements in a vertical plane deeper 
than ~ 50/mi in both geometries is due to the objective 
working distance that is about 200/mi. 

As a working fluid a buffer solvent [19] with sucrose 
concentration varied in the range between 47% (w/w)and 



67% (w/w) to tune viscosity in the range between rj s = 
0.012 Pa-s and rj s = 0.24 Pa-s at the working temperature 
of 22.5°C was used. 

To study the tumbling dynamics and statistics of poly- 
mer molecules in a flow 10 -3 ppm of A-DNA molecules, 
fluorescently labelled with YOYO-1 (Molecular Probes) 
at a dye/base ratio of 1:4 for ~ 1 hour, were added into 
the solvent. At equilibrium the coiled A-DNA has R g = 
0.73/mi, while the entire contour length is R ma x ~ 21/mi. 
So it may be considered as a "flexible" polymer with 
roughly 300 persistence lengths [1]. The maximum relax- 
ation time relevant for further analysis, r re i — 11±0.1 sec 
for a solvent with 62% (w/w) sucrose was found. r re i for 
solvents with different sucrose concentrations was taken 
proportional to its concentration variation [1,2]. Fluores- 
cently labelled DNA molecules were monitored via x63, 
1.4NA oil immersion objective (Zeiss) with 0.4 /zm depth 
of focus with a homemade inverted epi-fluorescent mi- 
croscope [4]. Images of the molecules were digitized, and 
their and <f> angles were automatically measured by ap- 
proximation of a molecule outline with an ellipse, which 
main axis was used as the end-to-end vector, R, to de- 
fine an inclination angle, either or cf> depending on the 
flow configuration (see Fig.lA,B,C). This approximation 
works well for sufficiently stretched molecules down to 
2 -T- 3/itm for the main axis. Two time series of different 
scenarios of tumbling that occur in a shear flow are pre- 
sented in Fig. ID. One mentions that tumbling can take 
a place via coiled as well as folded states. 

In the lower inset in Fig. 3 we present PDFs of for two 
values of Wi. Each plot is based on up to 30,000 points. 
It is clear seen that at low Wi = 1.6 PDF can be fitted 
by Gaussian almost entirely, while at high Wi — 17.6 
the PDF tails strongly deviate from the Gaussian. In 
Fig. 3 PDFs for three values of Wi are presented in the 
logarithmic coordinates. The PDF tails for all values of 
Wi decay algebraically with the exponent rather close 
to the theoretical, |#| -2 , indicated by a solid line on the 
plot. The larger Wi, the wider the scaling region at the 
tail as predicted by the theory [11]. Variations in Wi 
values were achieved by viscosity tuning as well as varia- 
tion in s by changing rotation speed between and 0.7 
sec -1 . The half-height width of PDFs, A9, decreases as 
oc Wi~°- 38±0 - 03 (see the upper inset in Fig. 3) in a rather 
good agreement with the recent numerical simulations 
based on the FENE model [14]. 

In the same set-up with the rotating rod we mea- 
sured PDFs of the tumbling time P(r), determined in 
the experiment as time between two following conse- 
quent conformations: two coiled, two folded, or coiled 
and folded ones (see Fig. ID). This tumbling time, r, is 
elongation based time defined by a certain length thresh- 
old R t h > R g , and it is different from that is angu- 
lar based. The lower the threshold, the higher r t [14]. 
We present in Fig. 4 P(t) for three values of Wi together 
with the exponential fits. Each plot is based on up to 600 
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points. The same data presented in the log- linear coor- 
dinates together with the fits demonstrate the fit quality 
(inset in Fig. 4). Figure 5 shows the characteristic tum- 
bling time, Tt, obtained from the slopes of the PDF tails, 
versus Wi. The value r t depends on the threshold value 
Rth/Rg- Wi was varied in two ways: full squares present 
the data, where Wi was tuned just by viscosity variation 
(it means the variation of T re ;), and open squares were 
obtained at different Wi by adjusting s. By presenting 
the same data as the ratio T re i/rt versus Wi in the inset 
in Fig. 5 we show, first, that r t is proportional to r re ; at 
lower Wi, and second, r re ;/r t grows as oc Wi°- 54 ± - 04 a t 
higher Wi in a rather close agreement with Eqs.(3,4). We 
would like also to point out the fact that the cross-over 
region for scaling in r t takes place at Wi ~ 10 in amaz- 
ing agreement between the experiment and the recent 
numerical simulations [14]. An obvious shift upwards of 
the data relatively to the simulations (inset in Fig. 5) is 
explained by about twice smaller value of Rth/Rg for de- 
termination of t in the simulations [14]. 

In Fig. 6 PDF of <f> measured in a micro-channel flow 
at Wi = 25 together with the fit by (sin0)~ 2 is shown. 
It is clear that the PDF maximum is located at (f> t ^ 
but limited angular resolution does not allow us to study 
its dependence on Wi. Similar data on P{4>), based on 
2000-3000 points for each plot, were obtained for several 
values of Wi. The half-height width of PDFs, A(j>, as 
a function of Wi is shown in the inset in Fig. 6 together 
with the fit and the results of the numerical simulations 
based on the FENE model [14]. The observed scaling 
oc Wi~°- 51±0 - 04 manifests the cross-over region between 
cx Wi~ x l z expected at high Wi and Wi~ x at lower Wi 
values (see Eq.(3,4)). 
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FIG. 1. A) Polymer orientation angles in a shear flow; B) 
Schematic diagram of rotating rod set-up; C) Schematic dia- 
gram of micro-channel set-up; D) Images of various polymer 
conformations consequent in time during tumbling. White 
bar is 3.5/im. 

FIG. 2. A) Azimuthal velocity vs vertical coordinate; in- 
set: the same vs radius. Region of measurements is around 
R = 700/im and at height of ~ 30^im. B) Squares are longitu- 
dinal velocity in horizontal plane, circles- the same in vertical 
plane. Region of measurements at 30 ± 10/im from the chan- 
nel wall and at height of ~ 50/im (dash line). Solid lines on 
all plots are theoretical. 
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FIG. 3. PDFs of \9\ in logarithmic coordinates: 
Wi = 1.6-pluses, Wi = 8.5-circles, Wi = 17.6-squares. A 
solid line is theoretically predicted asymptotic scaling for 
PDF tails. Lower inset: PDFs of 9 for Wi = 1.6-pluses, 
and Wi = 17.6-squares presented in log-linear coordinates; 
the solid curves are Gaussian fits. Upper inset: A9 vs Wi; 
squares-data, triangles-numerics. A solid line is the fit. 

FIG. 4. PDFs of r at several values of Wi. The solid lines 
are exponential fits. Inset: PDFs of r in log-linear coordinates 
with fits to the data at Wi = 4.7-triangles, Wi = 8.5-circles, 
Wi = 17.6-squares. 
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FIG. 5. r t vs Wi: full squares-data taken at constant s and 
varying viscosity; open squares-data taken at constant viscos- 
ity and varying s. The dash lines are guides to eye. Inset: 
T r& i/ T t vs Wi; triangles-numerics. Solid lines are the fits. 



FIG. 6. PDFs of (j> in log-linear coordinates at Wi = 25. A 
solid curve is the fit by oc (sin</>)~ 2 . Upper inset: Acj> vs Wi; 
squares-data, triangles-numerics. A solid line is the fit. 
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